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Abstract 

Molecular line data for the rotational transitions of a 
number of astrophysically interesting species have been 
collected and are made publically available through the 
www. These data, including energy levels, statistical weights 
Einstein A-coefficients and collisional rate coefficients, are 
essential input for non-LTE molecular radiative transfer 
programs. A computer program for performing statisti¬ 
cal equilibrium calculations is made publically available 
for use as well. This database should form an important 
tool in analyzing observations from current and future 
infrared and (sub)millimetre telescopes. 
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1. Introduction 

A wide variety of molecules has been detected in space 
to date, ranging from simple molecules like CO to more 
complex organic molecules including ethers and alcohols. 
Observations of molecular emission at millimetre and in¬ 
frared wavelengths, supplemented by careful and detailed 
modelling, are powerful tools to investigate the physi¬ 
cal and chemical conditions of astrophysical objects (e.g., 
Black 2000). For some objects, lines with a large range of 
critical densities and excitation temperatures are needed, 
since the densities can range from ~ 10 2 — 10 9 cm -3 and 
the temperatures from ~ 10 — 1000 K in the interstellar 
and circumstellar environments probed by current and 
future instrumentation. 

To extract astrophysical parameters, the excitation 
and radiative transfer of the lines needs to be calculated. 
A number of radiative transfer codes have been developed 
for the interpretation of molecular line emission [see van 
Zadelhoff et al. (2002) for a review]. The radiative transfer 
analysis requires molecular data in the form of energy lev¬ 
els, statistical weights and transition frequencies as well 


as the spontaneous emission probabilities and collisional 
rate coefficients. 

2. Molecular line data 

Literature data on the rotational transitions of 23 dif¬ 
ferent molecules are summarized and extrapolation of 
collisional rate coefficients to higher energy levels and 
temperatures are made. The data are made available for 
the community through the wwwQ Many of the data 
files presented here were adopted by Schoier et al. (2002) 
to model the circumstellar environment of the protostar 
IRAS 16293-2422. 

2.1. Energy levels and radiative rates 

The energy levels and Einstein A-coefficients are obtained 
from the JPL0 and CDMSQ catalogues. In some cases, 
better or more recent laboratory data have been used 
instead. Generally, only the ground vibrational state is 
retained, although in a few cases, e.g. CO and HCN, 
vibrationally excited levels are included. More study of 
collisional rate coefficients for vibrational transitions is 
urgently needed, as there are certainly cases where vi¬ 
brationally excited molecules are important. Examples 
include HC3N (Wyrowski et al. 1999) and torsionally ex¬ 
cited CH3OH. In the future, data files including excited 
vibrational states will be added for more molecules. 

Hyperfine splitting can be of importance in line trans¬ 
fer and introduce non-local effects for lines overlapping 
in frequency (e.g., Lindqvist et al. 2000). However, often 
the splitting between individual hyperfine components is 
small compared to the line-broadening, so that it can 
safely be neglected and treated as a single level for the 
purpose of excitation analysis. We present data files for 
some of the relevant molecules, such as HCN and CN, 
both with and without hyperfine splitting. 

1 http: //www. strw. leidenuniv.nl/~moldata 

2 http://spec.jpl.nasa.gov 

3 http://www.cdms.de 
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Figure 1. Calculated collisional de-excitation rate coeffi¬ 
cients for CO in collisions with para-H 2 for selected tran¬ 
sitions (symbols). Black lines indicate fits to the rate co¬ 
efficients of Flower & Launay (1985) (solid lines) and 
Flower (2001) (dashed lines). Also shown are the rate co¬ 
efficients presented by Schinke et al. (1985). 


2.2. Collisional rate coefficients 

The adopted collisional rate coefficients usually pose the 
largest source of uncertainty of the molecular data input 
to the radiative transfer analysis. Detailed summaries of 
the theoretical methods and the uncertainties involved 
in determining collisional rate coefficients are given by 
Green (1975) and Flower (1990). 

Published collisional rate coefficients cover only a lim¬ 
ited range of temperatures and energy levels, and extrap¬ 
olations are often necessary (see Fig. [I]) introducing ad¬ 
ditional uncertainties. Data files containing the originally 
calculated set of rates as well as files where an extrapo¬ 
lation has been carried out are available for all molecular 
species. 

In Fig. the predicted line intensities for the 15 low¬ 
est rotational transitions of CO of an isothermal homo¬ 
geneous sphere are shown using different sets of colli¬ 
sional rates. Deviations up to ~ 50% are found for the 
higher transitions which are sub-thermally excited. The 
lower transitions are close to LTE and less sensitive to 
the adopted set of collisional rate coefficients. 



Figure 2. Predicted CO line intensities using different sets 
of collisional rates for an isothermal homogeneous sphere 
with Tkin = 100 AT and tih 2 = 10 5 cm , compared to the 
values obtained using the CO-pH 2 rate coefficients from 
Flower (2001). Points connected with full lines are based 
on the original sets of rates while the dotted and dashed 
line styles indicate values obtained from fitting and ex¬ 
trapolation of the original rates 


3. Radiative transfer 

A radiative transfer code, RADE?^. is made available 
for public use as part of the data base. RADEX is a 
one-dimensional radiative transfer code aimed at solv¬ 
ing the statistical equilibrium equations using the escape 
probability formulation. The code makes no assumption 
for the geometry or large scale velocity fields and its 
basic assumption is that of an isothermal and homoge¬ 
neous medium. RADEX provides a useful tool in rapidly 
analysing a large set of observational data providing con¬ 
straints on physical conditions, such as density and ki¬ 
netic temperature. RADEX provides an alternative to the 
widely used rotation temperature diagram method (e.g., 
Blake et al. 1987; Goldsmith & Langer 1999) which re¬ 
lies upon the availability of many transitions of optically 
thin emission lines and is useful only in roughly constrain¬ 
ing the excitation temperature in addition to the column 
density. 

Benchmarking of the RADEX code is presented in 
Fig. I where various expressions for the escape probability 

4 http: //www. strw. leidenuniv.nl/~moldata/radex .html 
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Figure 3. Comparison of the predicted line strengths for 
the 10 lowest rotational transitions of CO for a homo¬ 
geneous isothermal sphere, with nn 2 = 10 5 cm -3 and 
Tkin = 50 K, using different methods. Upper panel - The 
total optical depth through the sphere at line centre r and 
excitation temperature T ex as a function of the upper ro¬ 
tational level J involved in the transition. Middle panel 
- The radiation temperature Tr obtained for each transi¬ 
tion using RADEX with different prescriptions of the es¬ 
cape probability, /3, and compared with the result from the 
Monte-Carlo code (MCC) of Schoier (2000). Also shown 
are the results for optically thin emission in LTE. Lower 
panel - Tr obtained from RADEX compared with the re¬ 
sults from MCC, <5Tr = T^ cc — Tr. 


are used. The output is compared both to an optically 
thin LTE analysis (rotation diagram method) and a full 
radiative transfer analysis using a sophisticated non-LTE 
code based on the Monte-Carlo method (MCC; Schoier 
2000). The molecular data used in this example are the 
collisional rate coefficients for CO collisions with para- 
112 from Flower (2001). Levels up to J = 29 are included. 
The excitation temperatures of the lines vary from nearly 
thermalized, for transitions involving low J-levels, to sub- 
thermally excited for the higher-lying lines. The optical 
depth in the lines is moderate (~ 1 — 2) to low. It is 
seen that the expressions of the escape probability for a 
sphere and the mean [(1 — e _T )/r] give almost identical 
solutions and close to that obtained from the full radiative 
transfer. However, the slab geometry gives quite different 
results, in particular for high lying lines. The optically 
thin approximation, where the gas is assumed to be in 
LTE at 50 K produces the largest discrepancy and only 
gives the correct answer for the J = 1 —» 0 line, where 
the requirements are fulfilled. 
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